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ABSTRACT 

We report the discovery of the likely white dwarf companions to radio millisecond 
pulsars 47Tuc Q and 47Tuc S in the globular cluster 47Tucanae. These blue stars 
were found in near-ultraviolet images from the Hubble Space Telescope for which we 
derived accurate absolute astrometry, and are located at positions consistent with 
the radio coordinates to within 07016 (0.2 a). We present near-ultraviolet and optical 
colors for the previously identified companion to millisecond pulsar 47 Tuc U, and we 
unambiguously confirm the tentative prior identifications of the optical counterparts 
to 47 Tuc T and 47 Tuc Y. For the latter we present its radio timing solution for the 
first time. We find that all five near-ultraviolet counterparts have [/300 — ^390 colors 
that are consistent with He white dwarf cooling models for masses ^ 0.16 — 0.3 Mq 
and cooling ages between ^0.1 — 6 Gyr. The Ha — 7^625 colors of 47 Tuc U and 47 Tuc 
T indicate the presence of a strong Ha absorption line, as expected for white dwarfs 
with a H envelope. 

Key words: Binaries: general - globular clusters: individual (47Tucanae) - pul¬ 
sars: general - pulsars: individual (PSR J0024-7204Q, PSR J0024-7204S, PSR J0024- 
7204T, PSR J0024-7203U, PSR J0024-7204Y) 


1 INTRODUCTION 

Millisecond pulsars (MSPs) are rapidly spinning neutron 
stars with spin periods of 20 ms or less, magnetic field 
strengths B ^ 10®“^ G, and charact eristic spin-dow n ages 
Tc ~ 10^“^° yr (see e.g. the review bv iRanso r^ l2008lV Their 
proper ties can be explain ed by the recycling formation sce¬ 
nario (jAloar et al.l Il982l ) , in which an old neutron star is 

* E-mail: l.e.riverasandoval@uva.nl 

^ Based on proprietary and archival observations with the 
NASA/ESA Hubble Space Telescope, obtained at the Space Tele¬ 
scope Science Institute, which is operated by AURA, Inc., under 
NASA contract NAS 5-26555. 


spun up by the accretion of mass from a companion star in 
a low-mass X-ray binary (LMXB) phase. 

Up to now, more than 200 MSPs have been discovered 
in the Galaxy, and a large fraction of these reside in globular 
clusters (GCs). The number of MSPs per unit mass is higher 
in GCs tha n in the Gal actic disk by two to three orders of 
magnitude iln eirell201^ . In the context of the recycling sce¬ 
nario, this overabundance is no surprise. The LMXB pro¬ 
genitors of MSPs have long been known to occur in GCs 
with a much higher frequency than in the rest of the Galaxy 
Ja a [mi) as a result of the high stellar encounter rates 
that favor their formation in dense clusters. 

Most MSPs are found in binary systems. Studies of 
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Table 1. Parameters of the millisecond pulsars in 47Tuc whose companions are discussed in this paper, for a reference epoch MJD 
51600. Positions, spin periods (Ps), orbital periods (Pfo), limits on companion masses (me; for an assumed neutron star mass of 1.35 
Mct)) . and characteris t ic a^e s (tc) are from Freire et al. (in preparation). X-ray luminosities in the 0.1 — 10 keV band (Lx) were taken 
from lBogdanov et al.l (l2006lV Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, areminutes, 
and arcseconds. Numbers in parentheses are Icr uncertainties in the last significant digit, errors in Ps are in thirteenth decimal, while 
errors in are in the ninth decimal or smaller. Values of Tc are lower limits based on the 2 a upper limit for the spin-period derivatives. 


MSP 

R.A. (J2000) 

Decl. (J2000) 

Ps 

(ms) 

Pb 

(days) 

ruc 

(Me) 

Tc 

(Gyr) 

Lx 

(lO^o erg s-i) 

Q 

00:24:16.4903(2) 

-72:04:25.1653(8) 

4.03 

1.189 

> 0.18 

> 1.43 

3.9 

S 

00:24:03.97945) 

-72:04:42.35355) 

2.83 

1.201 

> 0.09 

> 0.91 

8.9 

T 

00:24:08.5487(7) 

-72:04:38.928(3) 

7.58 

1.126 

> 0.17 

> 0.32 

2.9 

U 

00:24:09.83626(7) 

-72:03:59.6889(3) 

4.34 

0.429 

> 0.12 

> 1.91 

5.1 

Y 

00:24:01.4016(3) 

-72:04:41.837(1) 

2.19 

0.521 

> 0.14 

> 2.2 

3.9 


the orbital parameters and the companion properties can 
provide valuable information on the ad vanced evolution of 
LMX Bs, the pulsar spin-up process (e.g. ISmarr BlandfordI 
Il976l ). and (for MSPs in GCs) the effects of dynamical en- 
cou nters on the characteris tics of the MSP population (see 
e.g. lVerbunt fc Freirell201 j ). Various classes of binary MSPs 
can be distinguished; each of these are found in the Galac¬ 
tic field and in GGs. The systems that are expected to re¬ 
sult from the standard recycling scenario have low-mass (a 
few tenths of M©) white dwarf (WD) companions that de¬ 
scended from the original mass donors in the LMXBs. Most 
MSPs are found in this configuration. The very-low mass bi¬ 
nary MSPs (or black widows) likely represent the very latest 
stages of evolution, in which the companion (with only a few 
hundredths of M© in mass remaining) is being ablated by 
the pulsar. Finally, the so-called redback MSPs are charac¬ 
terized by radio ec lipses and non -degenerate companions of 
0.2 — 0.4 M© (see lRobertsll2013l for a review). These may 
be systems in transition between an LMXB phase and a 
detached phase where the MSP can be detected as a ra¬ 
dio pulsar, or (in globular clusters) the result of exchange 
encounters. 

The identification of MSPs at other wavelengths en¬ 
ables a range of studies into the physical processes that 
govern the evolution of the MSP binaries. In GGs, the spa¬ 
tial resolution of the Chandra X-ray Observatory and the 
Hubble Spaee Teleeope (HST) is a necessity to overcome the 
crowding. Indeed, the first unambiguous X-ray identifica¬ 
tion of MSPs in globular clusters was made with the first 
Chandra observation of 47Tuc, in which fifteen of the 23 
MSPs the n known were first de tected and precisely located 
in X-rays (|Grindlav et akllioOlh . So far, HST counterparts 
to globular-cluster MSPs have been found for only twelve 
syste ms, two of which are ten tative identifications (47 Tuc T 
and Y: lEdmonds et al.ll2003a). The other ten sys tems in GGs 
are 47 Tuc U and W (lEdmonds et al.l[200l|:l20Q2l). NGC 6397 


A (lEerraro et al.l 200ll ), NG G 6752 A (Bassa et al.] 120031 . 

lEerraro et al.l 20031). M4 A (|Sigurdsson e^ . l2003|). NGG 


6266 B (iGocozza et al.l 20081). M28 H and I ( Pallanca et ^ 

I2OIOI. I2OI3I), M5 G |Pallanca et Zl l2014l ) and M71 A 
(|Gadelano et al.l l2015l ). Glearly, more identifications are 
needed to get a better understanding of the possible dif¬ 
ferences between MSP properties in GGs and in the field. 

Right after Terzanb, 47Tucanae (47 Tuc) is the GG 
where most radio MSPs ha ve been identified so far : a to- 
tal of 23 were discovered ([Manchester et al.] Il990l . Il99ll : 


iRobinson et al.l [l995l : JCamilo et ^l2000l ). fifteen of which 
are in binary systemqj. To date, accurate radio-timing po¬ 
sitions have been derive d for thirteen of the 47 Tuc binary 
MSPs ([Ereire et al.ll20o3 . and in preparation), which allows a 
search for their counterparts at ultraviolet and optical wave¬ 
lengths. So far, secure identifications have been found for 
only two of them. The companion of PSR J0024-7204W, 
or 47 Tuc W in short, is a non-deg enerate star that may 
be (close to) filling i ts Roche lobe ([Edmonds et al.l I 2 OO 2 I : 
iBogdanov et ahlboOSl ): it was among the first MSP binaries 
to be classified as a ’’redback” system. The companion of 
47 Tuc U, on the other hand, is typical for a system that is 
the outcome of the canonical recycling scenario. Edmonds 
et al.(2001; EOl hereafter) found that it is likely a He WeEI 
with a mass of ~0.17 . lEdmonds et al.l (2003a; from now 

on E03a) suggested identifications for two more systems, 
47 Tuc 43 and T, but from their data they could only de¬ 
rive upper limits on the U — V colors of the proposed blue 
counterparts. 

We have carried out a search for counterparts to the 
47 Tuc MSPs using new near-ultraviolet (NUV) imaging 
data obtained with the HST Wide Eield Gamera 3 (WEG3), 
and optical data (including deep Ha imaging) from the Ad¬ 
vanced Gamera for Surveys (AGS). The NUV data are par¬ 
ticularly useful for looking for WD companions, whose blue 
spectral energy distribution yields a high contrast against 
the typically red cluster stars. Here we present the results of 
our search, which includes the discovery of the counterparts 
to 47 Tuc Q and S, and the corroboration of the suggested 
identifications of 47 Tuc T and Y. We find that, like in 47 Tuc 
U, the MSP companion in all four systems is likely a low- 
mass (0.16 — 0.3 M©) He WD. 


^ http://www.naic.edu/~pfreire/GCpsr.html 
^ He WDs are low-mass (< 0.45 M©) white dwarfs with He-rich 
cores. They are thought to have formed as a result of a mass- 
transfer episode in a close binary before He ignition, removing 
much of the He WD-progenitor’s envelope and exposing its He- 
rich core. 

The optical counterpart to 47 Tuc Y was ori ginally reported as 
the counterpart to the Chandra source W 82 bv lE03a| . This source 
was classified as a candidate cataclysmic variable, as the radio 
position for 47 Tuc Y was unknown at the time. After the position 
became available, Bogdanov et al. 2006 reported its association 
with W 82, implying the detection of the MSP companion by 
E03a. 
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This paper is structured as follows. In Section 2 we de¬ 
scribe the HST observations, and the X-ray data used in 
this work. The HST data analysis is described in Section 3. 
In Section 4 we present the search for counterparts to the 
47Tuc MSPs in the HST data, and we show our results for 
the new identifications and the known WD companion of 
47Tuc U. Finally we discuss our results in Section 5. 


2 OBSERVATIONS 

The NUV data used in this worlfl belong to HST program 
GO 12950 taken with the UVIS channel of the WFC3. The 
data set consist of eight orbits divided in two consecutive 
visits on 2013 August 13 starting at 02:10 UTC and finish¬ 
ing at 11:01 UTC of the same day. Each visit started with 
one orbit of four F390W exposures of about 580 s each. This 
was followed by three orbits of twelve F300X exposures in 
total, with individual exposure times of about 609 s, result¬ 
ing in a total exposure time of 3881 s in the F390W filter 
and 14256 s in the F300X filter. The native scale of the 
WFC3/UVIS images is O'.'04 pixeP^ The UVIS detector is 
composed of 2 chips, having a chip gap between them equiv¬ 
alent to 172 on the sky. In order to improve our spatial res¬ 
olution, we used a 4—point dithering pattern that includes 
fractional-pixel offsets. In both visits an offset between im¬ 
ages was imposed to ensure coverage of the chip gap. The 
im ages were centered nea r the center of 47Tuc as reported 
bv iGoldsburv et ^ (|201lh . Given that the WFC3/UVIS has 
a field of view (FOV) of 162" x 162", our images only cover 
~ 23% of t he area insid e the half-mass radius of the cluster 
{ih = 3.'17, [HarUi[l9^, 2010 editior0) 

The optical observations reported here were taken with 
the Wide Field Channel (WFC) on the ACS under pro¬ 
gram GO 9281. We obtained images in three filters: F435W, 
F658N, and F625W, which encompasses Ha and the sur¬ 
rounding continuum. The observations were spread over 
three visits on 2002 September 30 (01:25-06:31 UTC), Oc¬ 
tober 2/3 (20:37-00:51 UTC), and October 11 (02:04-10:40 
UTC). The visits were scheduled to be simultaneous with 
the deep 47 Tuc Ch andra observations that are presented in 
iHeinke et al.l (|2005l : hereafter H05). The total exposure time 
amounted to 955 s in F435W, 1320 s in F625W, and 7440 s 
in F685N. Observations through all three filters were taken 
in each visit, resulting in near-simultaneous colors. Since X- 
ray emitting binaries frequently display variability in X-ray 
and optical wavelengths, this provides an optimal estimate 
of the system’s actual colors, which is unaffected by long¬ 
term variability in one or both wavelength bands. The FOV 
of the WFC is 202" x 202", imaged with two CCD detectors 
separated by a 2745-wide gap. We used a dither pattern 
with fractional-pixel offsets to be able to improve the native 
WFC pixel scale of 0705 pixel” ^ during the image stacking. 
A complete description of the GO 9281 data and results is 
forthcoming (van den Berg et al. in preparation). Of the 


^ We have not analyzed the existing STIS far-ultraviolet data of 
the center of 47 Tuc, because the five objects studied in this paper 
are not included in the small field of view of that data set. 

® Catalog of parameters for Milky Way GCs, 
http: //www.physics.mcmaster.ca/Fac_Harris/mwgc.dat 


twenty MSPs in 47 Tuc with radio-timing positions, nine¬ 
teen are covered by the UVIS and WFC images; only 47 Tuc 
X, a binary MSP, was missed. 

In this work we make use of the catalog of X-ray sources 
in 47 Tuc from H05 (see that paper for details about the data 
reduction), which was obtained from deep observations (281 
ks in total) with Chandra. All nineteen MSPs covered by our 
NUV/optical data are matched to an X-ray source in this 
catalog. 


3 DATA REDUCTION AND ANALYSIS 

The data reduction of our UVIS data starts with the flat- 
fielded images produced by the standard UVIS pipeline 
(CALWF3 version 3.1.2). To correct for the charge transfer 
efficiency (CTE) degradation of the UVIS detectors, we used 
the CTE correction softwar^H provided by the Space Tele¬ 
scope Science Institute (STScI). The resulting images still 
suffer from geometric distortion, which we correc ted for with 
the DrizzlePac0 software (Gonzaga et al. l2012h . Astromet¬ 
ric alignment of the individual images was performed with 
the TWEAKREG task in Drizzlepac. Subsequently, the AS- 
TRODRIZZLE task was used to combine the aligned images 
into one stacked master frame in each filter that is distortion 
free, and cleaned of cosmic rays and detector artifacts like 
hot pixels or bad columns. The resulting master images in 
E390W and E300X are twice-oversampled to a pixel scale of 
0702 pixel”^ 

The reduction steps for the WEC images are similar to 
those for the UVIS data, except that the images produced 
by the ACS calibration pipeline (CALACS version 8.2.0) 
are already corrected for CTE losses. The stacked master 
images in E435W, E658N, and E625W have a final pixel 
scale of 07025 pixel”^. 


3.1 Astrometry 

The radio coordinates of the 47 Tuc MSPs were derived us¬ 
ing the Solar System ephemeris DE/LE 405. This is speci¬ 
fied in the International Celestial Reference Erame (ICRE), 
which is aligned to the International Celestial Reference Sys¬ 
tem (ICRS). To facilitate the identification of NUV/optical 
counterparts we also aligned our UVIS and WEC images to 
the I CRS using stars in the UCAC2 catalog (jZacharias et al.l 
l2004h . Since isolated and unsaturated UCAC2 stars are 
scarce in the HST images of the cluster core, we derived 
secondary astrometric standards from a ground-based im¬ 
age of 47 Tuc as an intermediate step. We obtained from 
the ESC archive a 30-s V image taken with the 2.2m/Wide 
Eield Imager (WEI) at La Silla, Chile on 2002 October 29. 
An astrometric solution of the part of the image that con¬ 
tains the cluster center was derived using 225 UCAC2 stars; 
fitting for zero-point offset, rotation, pixel scale, and distor¬ 
tions gives rms residuals of about 07035 in right ascension 
and declination. We then selected 126 secondary standards 
from the WEI image to calibrate a 10 s E435W exposure 


http://www.stsci.edu/hst/wfc3/tools/cte_tools 
® http://drizzlepac.stsci.edu 
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Table 2. Astrometric and photometric results for the MSP companions discussed in this paper. Errors in the absolute astrometry of 
the celestial positions are 0''074 (Icr). Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, 
arcminutes, and arcseconds. Columns 4 to 8 give the calibrated photometry on the Vega-mag system and the errors from DAOPHOT. 
The project ed distance of the com panions from the cluster center ((5c) was calculated using the coordinates of the center of 47Tuc as 
reported bv IColdsburv et al.l (|201l|j . The offsets between the X-ray and NUV positions (last column) were determined using the positions 
from H05. The value of ax is the combined error in the NUV and X-ray positions. 


MSP 

companion 

R. A. 
(J2000) 

Decl. 

(J2000) 

U300 

-P390 

-P435 

Rq 25 

Ha 

(arcsec) 

Radio 

offset 

(arcsec/cr) 

X-ray^ 

source 

X-ray 

offset 

(arcsec/(Tx) 

Quv 

00:24:16.493 

-72:04:25.15 

22.95(3) 

23.41(4) 

23.84(6) 

23.6(1) 

22.93(8) 

56.89 

0.012/0.2 

W104 

0.05/0.6 

Suv 

00:24:03.976 

-72:04:42.34 

23.16(3) 

23.77(5) 

23.47(5) 

- 

- 

13.09 

0.016/0.2 

W77 

0 . 342 / 5.7 

Tuv 

00:24:08.551 

-72:04:38.92 

23.04(2) 

23.72(3) 

23.9(1) 

23.5(2) 

23.4(1) 

19.02 

0.011/0.1 

W105 

0.25/2.5 

Uuv 

00:24:09.836 

-72:03:59.69 

20.545(9) 

20.98(1) 

20.860(5) 

20.82(1) 

20.84(1) 

56.33 

0.002^.03 

Wll 

0.07/0.9 

Yuv 

00:24:01.402 

-72:04:41.84 

22.09(2) 

22.51(5) 

22.36(5) 

- 

- 

22.66 

0.002/0.03 

W82 

0.21/1.9 


^ Using the nomenclature in H05. 

^ The X-ray position used for 47 Tuc S corresponds to the position of a single detected X-ray source that is the combination of X-ray 
emission from 47Tuc F and 47 Tuc S. This explains the apparent large X-ray astrometric offset. 


from our GO 9281 program, which resulted in rms residu¬ 
als of about 0^.^027 in each direction. Finally, this solution 
was transferred to our stacked WFC images with negligible 
errors using >10 000 stars, and to our stacked UVIS im¬ 
ages with rms residuals of about 0V027 in each direction 
using >11700 stars. The final 1 sigma error in the absolute 
astrometry of our images is the quadratic sum of all the 
aforementioned errors, plus a term that represents the error 
with which the UCAC2 is tied to the ICRS. The latter is 
dominated by the estimated systematic error in UCAC2 po¬ 
sitions of 10 mas (jZacharias et al.|[2004l ). For the WFC and 
UVIS images this results in an error of 0V064 and 0V074, 
respectively. Errors in the radio positions are negligible. 

H05 aligned the astrometry of their X-ray source cat¬ 
alog to the radio positions of seventeen MSPs detected by 
Chandra. Therefore, we can indirectly check the alignment 
between our NUV and radio positions by comparing the 
NUV and Chandra astrometry. We have done this by com¬ 
paring the NUV and X-ray positions of se venteen sources 
with published finding charts in E03a and lEdmonds et al.l 
(2003b). We find that the NUV/X-ray offset is very small, 
about -0V023zb0V008 in right ascension and -0V003zb0V009 
in declination in both filters, so about one oversampled pixel. 

3.2 Photometry 

Given that the images of the core of 47 Tuc are very crowded, 
we carried out point spread function (PSF) photometry for 
the stacked UVIS and WFC frames by using standard IRAF 
and DAOPHOT tasks. We used isolated stars to build a 
variable PSF for each filter. The obtained PSF was fitted 
to all detected stars in each stacked frame. Photometric cal¬ 
ibration of the extracted magnitudes was performed onto 
the Vega-mag system using the zero points given at STScFs 
web page^; for UVIS we used the zero points provided for 
0^.^4 apertures. In the remainder we denote calibrated mag¬ 
nitudes in the F300X, F390W, F435W, F625W, and F658N 
filters by U300, H390, R435, ^625, and Ra, respectively. The 
absolute magnitudes were calculated for an adopted distance 
modulus (m — M)o = 13.36 zb 0.02 zb 0.06, corresponding to 

® http://www.stsci.edu/hst/wfc3/phot_zp_lbn 
http: / / WWW. stsci. edu/hst / acs / analysis / zeropoints 


the 4 7 Tuc distance of 4.69 zb 0.04 zb 0.13 kpc (|Woodlev et al.l 
l2012h . The first error term corresponds to the random er¬ 
ror, and the second to the systematic error of the distance 
determination. 

The UVIS data were also independently analyzed using 
software based on the program developed fo r the ACS Cloh- 
ular C luster Treasury Project, described in lAnderson et al.l 
(|2008h . which gave consistent results with the DAOPHOT 
analysis. The reductions were done using the latest version of 
this software known as KS2. Star finding required that stars 
be detected in both the F300X and F390W frames. Pho¬ 
tometry was performed by PSF fitting using library PSFs 
with perturbations to better match the specific PSF in each 
exposure. The final KS2 photometry comes from averaging 
the fluxes from the individual images, with sigma clipping 
applied to remove outlying values due to cosmic rays, de¬ 
fective pixels, etc. A total of 148717 stars were detected in 
the region covered by the UVIS imaging. Color-magnitude 
diagrams were constructed to choose the best KS2 options 
for photometry of the faint stars (e.g. the fitting radius), 
from the tightness of the fiducial sequences, and to iden¬ 
tify the counterparts for the MSPs. Drizzle-combined mo¬ 
saic images of this region were also produced using the STScI 
PyRaF routine ASTRODRIZZLE from the Drizzlepac pack¬ 
age. The drizzle-combined images were oversampled by a 
factor of two, in order to increase the effective resolution. 
Calibration of the KS2 photometry to the Vega-mag system 
was performed by doing aperture photometry on moderately 
bright, isolated stars in the image mosaics within a C!l ra¬ 
dius aperture, finding the aperture correction to an infinite 
radius aperture from ISR WFC3 2009-3^^. calculating the 
median offset between the KS2 photometry and the aperture 
photometry for these stars, and applying the zeropoints from 
the i/S'T WFC3 zero point website®. 

To take into account the effect of extinction on the UVIS 
photometry, we converted t he reddening towar ds 47 Tuc of 
E{B -V) = 0.04 ± 0.02 llSalaris et alJ l2007tl to the ex- 
tinction in the F300X and F390W filters using the UVIS 
Exposure Time Calculator of HST. We obtained extinction 
values of A(F300X)=0.26zb0.13 and A(F390W)=0.18zb0.09, 


http://www.stsci.edu/hst/wfc3/documents/ISRs/WFC3- 
2009-38.pdf 
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respectively. To compute the e xtinction in the WF C fil¬ 
ters, we used the conversions in ISirianni et al.l (|2005l i. This 
yields A(F435W)=0.16±0.08, A(F658N)=0.10±0.05, and 
A(F625W)=0.11±0.05. 


4 RESULTS 

The search for MSP counterparts was done in the astro- 
metrically calibrated NUV images. These frames are less 
crowded than the optical images, and contain fewer bright 
and saturated stars that could affect the detection of faint 
counterparts. We looked for counterparts within 3 a (0^22) 
of the radio positions. 

We found blue stars inside the match circles of the bi¬ 
nary MSPs 47Tuc Q, S, T, U, and Y. These objects are 
excellent astrometric matches with the radio positions, with 
offsets that are smaller than 0^016 (0.2 a). All other stars in¬ 
side the 3 a error circles lie further from the radio positions, 
and do not have blue Usoo — B 390 colors. 

The blue stars that match with 47Tuc T, Y and U 
are the likely WD coun terparts proposed by EOl, E03a and 
iBogdanov et al ] (l2006ll . while for 47Tuc Q and S these are 
new discoveries. Based on the similar locations in the B 390 
versus f /300 — B 390 color magnitude diagram (CMD), and 
the good alignment with the radio positions, we consider the 
stars near 47 Tuc Q and S to be the likely counterparts, and 
WD companions, to these MSPs, as well. A more detailed 
discussion of the individual systems is given in Sections 14.11 
to 14.51 In the remainder of the text we use capital letters 
(without subscript) to refer to the MSPs detected at radio 
wavelengths, while we add the subscript ”UV” to refer to the 
NUV counterparts, e.g. 47 Tuc Q versus Quv- Binding charts 
in the E300X and E390W filters are presented in Eigure [T] 
Eigures [2] and [3] are NUV and optical CMDs for 47 Tuc that 
show the photometry for Quv, Suv, Tuv, Uuv and Yuv- 
In Table [1] we summarize the properties of the five MSPs 
as determined from radio and X-ray observations, while our 
astrometric and photometric results can be found in Table 

Q] 

We have estimated the parameters of Quv, Suv, Tuv, 
Uuv and Yuv assuming that they are He WDs. In Eigures [2] 
and [ 3 ] we have included new evo lutionary tracks for H e WD 
cooling models as described in ISerenelli et al.l (|2Q02l b but 
with a metallicity as appropriate for 47 Tuc of Z = 0.003 
{[Fe/H] = —0.78. iThvgesen et al .11 20 i'l l. Masses were esti¬ 
mated by determining the mass range that corresponds to 
the two cooling tracks that bracket the location of the com¬ 
panion in the NUV CMD, including 1 cr error bars. Errors 
include the DAOPHOT uncertainties, errors in the calibra¬ 
tion, and in the adopted distance modulus and extinction 
(see Section [3.21) . The corresponding ranges in age, effective 
temperature Te//, surface gravity log^, and bolometric lu¬ 
minosity L were determined by interpolating linearly along 
the two bracketing tracks using the absolute magnitude in 
the E390W filter (Mb^qo) as the independent variable. Table 
[3] summarizes the results. 

Cooling ages derived from our models (column 3 in Ta¬ 
ble [3|) represent the time since the point of highest temper¬ 
ature in the cooling track until the observed location along 
the track. Prior to this time lapse is the time between Roche- 
lobe detachment (when the mass transfer from the He WD 



Figure 1. Finding charts of the companions to MSPs 47 Tuc Q, 
S, T, U and Y. Images in the filter F300X are shown on the left 
side while images in F390W are shown on the right. The solid red 
lines show the 3 a match circles centered on radio positions of the 
pulsars, which take into account the astrometric UVIS errors. The 
dashed cyan circles are the 3 a match circles of Chandra sources. 
The companion candidates are indicated by solid orange circles. 
Each image is I70 x I70 in size. 
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- 2-10 1 2 

^U300~ ^B390 


Figure 2. C /300 — ^390 versus S 390 CMD based on DAOPHOT photometry of our UVIS images. Each MSP companion is denoted with 
a letter and a filled magenta circle. The green filled circle represents the blue object detected by DAOPHOT inside the 3 cr match circle 
of 47Tuc I. The orange filled circle is the same blue object detected by KS2 insid e the 3 cr match circle of 47Tuc I. Calibrated magnitudes 
were converted to absolute magnitudes using a distance to 47Tuc of 4.69 kpc (| Woodley et al.ll2012|j and reddening E{B — V) = 0.04 
JSalaris et ahll^OTlj . The lines correspond to cooling models for He WDs with thick H envelopes as in ISerenelli et al ] ll 20 n 2 l) . computed 
for an initial metallicity of Z = 0.003. Labels at the top of each track represent the mass of the He WD in units of solar mass. 


progenitor to the neutron star stops) until the point where 
the highest Te// is reached. This phase lasts only a few Myr 
for masses 0.25 Mq and higher, but can amount to ~ 1.5 Gyr 
for the 0.16 Mq model, ~1 Gyr for 0.175 M©, and ~500 Myr 
for 0.2 Mq. The exact value depends on the uncertain de¬ 
tails of t he binary evo l ution , and sets a lower limit to the 
total age. Ilstrate et al.l (|2014|j have derived a general expres¬ 


sion for the duration of this ”proto-WD” phase (Atproto) in 
terms of the He WD mass, which for our MSP companions is 
given in the last column of Table [S] We find that the proto- 
WD time scale in our models and the value of Atproto are in 
general agreement. 

We have estimated the number of He WD-like objects 
that we expect to find by chance within a separation of 
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Figure 3. S435 — i ?625 versus Rq 25 (left) and Ha — Rq 25 versus i ^625 (right) CMDs from our WFC DAOPHOT photometry. MSP 
counterparts are plotted with filled re d circles. Calibrated Vega-mag magnitudes were conver ted to dereddened absolute magnitudes 
using a distance to 47Tuc of 4.69 kpc (|Woodlev et al .1120121 ') and reddening E{B — V) = 0.04 (|Salaris et al.ll20o'^ h We include He WD 
cooling tracks (Serenelli et al. 2002) for masses in the range 0.16 - 0.45 Mq and an initial metallicity of Z = 0.003. Photometry errors 
include the DAOPHOT uncertainties, and errors in the distance modulus and E(B — V). 


0('016 (i.e. the largest offset among the suggested matches 
for 47Tuc Q, S, T, U, and Y; see Table [2]) from the radio 
positions. To this end, we determined the density of blue 
objects (stars arcsec“^) enclosed by the cooling tracks for 
masses 0.16 ^ MwdIMq ^0.30 and absolute magnitudes 
6 ^ Mbsqo ^11 (see Figure [2]). For the total area of the 
3 a match circles of nineteen MSPs, the expected number is 
very small, 9.5 x 10“^. This suggests that our five identified 
MSP counterparts are indeed not chance coincidences. 

For the binary MSP 47Tuc I we find a very blue ob¬ 
ject at 0('12 (1.6 cr) with a U300 — B390 color that puts it 
on or very close to the CO WD sequence (Figures [2] and 2]). 
This offset is much larger, and its color much bluer, than 
those for the five binary MSPs discussed above. We esti¬ 
mate the probability that this source is a spurious match 
by considering that blue {U300 — B390 < —0.65) objects are 
found in the error circles of two (47Tuc D and N) of the 
seven isolated MSPs in the UVIS field. Due to their intrinsic 
faintness in the NUV/optical, no counterparts are expected 
for the isolated MSPs. The blue objects close to 47Tuc D 
(at 0Vl4 or 1.8 cr, and at (-1.44, 11.50) in the CMD of Fig¬ 
ure [2]) and 47Tuc N (at 0 VI 6 or 2.1 cr, and at (-0.65, 9.68) 
in the CMD)—are therefore almost certainly chance align¬ 
ments. The probability to find a blue star in the 3cr error 


circle by chance in a single trial is therefore, approximately, 
2/7 ~ 0.29. This results in a (binomial) probability to find 
a blue random match for exactly one of ter0 binary MSPs 
in our field of ~14%. While we consider the NUV identifica¬ 
tion of 47 Tuc I not very likely, we give more details on this 
object in Section 4.6, and provide a finding chart in Figure 

m 

For 47Tuc E, H, J (binaries), and L (isolated) not a sin¬ 
gle NUV object was detected inside their 3 a match circles. 
For 47 Tuc O, R, and W (binaries) and C, F, G, and M (iso¬ 
lated) the objects in the match circles have colors that place 
them on the main sequence or sub-giant branch of 47 Tuc 
and lie at least OVll (1.4 cr) away from the radio positions. 
Without any additional information we cannot conclusively 
say if any of these are likely counterparts, and we do not 
discuss them further in this paper. We not e that the true, 
non-d egenerate, counterpart to 47 Tuc W ([Edmonds et al l 
I2OO2I ) is not detected in the NUV images. This is not un¬ 
expected as it is relatively faint (mean V = 22.3) and red 
(compared to Quv, Suv, Tuv, Uuv, and Yuv). 


11 For 47 Tuc U and W excellent counterparts have already been 
identified. 
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Table 3. Limits on companion masses {rric] for an assumed neutron star mass of 1.35 Mq) as derived from the mass function (Freire 
et al. in preparation). Masses (Mwd), cooling ages, effective temperatures (Tg^^), surface gravities (log g) and luminosities (L) for the 
MSP companions studied, as derived from our NUV photometry and our Z = 0.003 He WD cooling models. In column 8, the inclination 
angle (i) as inferred from the mass function and the estimated masses Mwd. The last column gives At proto, which is the lapse of 
time from th e Roc he-lobe detachment until the proto-He WD reaches the WD cooling track, calculated using the relation derived by 
list rate et al.l (|2014l L 


Companion 

ruc (Mq) 

Mwd (Mq) 

Age^ (Gyr) 

Te// (kK) 

log 9 

log L (Lq) 

i (°) 

Atproto 

(Gyr) 

Quv 

> 0.18 

0 . 175 - 0.20 

5 . 2 - 0.3 

8.9 - 10.6 

6 . 4 - 6.7 

- 1.9 - - 1.8 

90-63 

1 - 

0.4 

Suv 

> 0.09 

0.20 - 0.25 

0 . 4 - 0.1 

9 . 7 - 11.3 

6 . 7-7 

- 2.0 - - 1.9 

28-22 

0 . 4 - 

- 0.1 

Tuv 

> 0.17 

0.20 - 0.30 

0 . 4 - 0.1 

9.9 - 12.8 

6 . 7 - 7.2 

- 2.0 - - 1.8 

58-36 

0 . 4 - 

0.02 

Uuv 

> 0.12 

0.16 - 0.1752 

0 . 2 - 0.3 

10.1 - 12 

5 . 5 - 5.8 

- 0.9 - - 0.8 

51-46 

1 . 9 - 

- 1.3 

Yuv 

> 0.14 

0.16 - 0.175 

2.1 - 1.8 

8.9 - 10 

6 - 6.2 

- 1.6 - - 1.5 

60-53 

1 . 9 - 

- 1.3 


^ Ages represent the time between the point of highest Tg// on the cooling track until the observed point along the track. 

^ An upper mass limit of 0.17 Mq was found by EOl, for which they estimated an age of ~ 0.6 Gyr, T^ff = 11 000 K and log g = 5.6 . 


4.1 47Tuc Q 

In the CMD of Figure O the counterpart of 47Tuc Q is lo¬ 
cated between the main sequence of the Small Magellanic 
Cloud (SMC) and the CO WD cooling sequence, at a sim¬ 
ilar 7/300 — ^390 color as Uuv, but about about 2.4 mag 
fainter in B 390 . E03a reported the presence of a variable 
main-sequence turnoff star very close (at 0('24) to the posi¬ 
tion of 47Tuc Q. This star, which they labelled ”nQ”, is the 
bright object to the south-west of Quv in Figure [T] They 
discarded nQ as the true counterpart to the MSP, since its 
period of optical variability does not match the orbital pe¬ 
riod of the pulsar derived from radio timing. The faint blue 
star Quv lies closer to the radio position than nQ (0('012 or 
0.2cr), but was not detected by E03a. It lies in the PSF wings 
of two stars (nQ and a star on the 47 Tuc sub-giant branch), 
which makes Quv difficult to detect at optical wavelengths. 

In our Re25 versus B435 — ^625 CMD, Quv also lies 
in the region where He WDs are expected. But unlike the 
counterparts to 47 Tuc U and T that lie to the red side (Ha 
faint side) of the main sequence in the Rq 25 versus Ha — Rq 25 
CMD (see Sect. 14.31 and 14.51 and Figure [3|), there is no in¬ 
dication for the presence of a broad Ha absorption line in 
Quv: its Ha —7^625 color is consistent with the color of most 
other stars in the field with similar R 625 magnitude. If Quv 
actually is a He WD, the Ha —7^625 color could be explained 
if there is some residual Ha emission in the system that ef¬ 
fectively fills in the Ha absorption line of the H-rich outer 
layer of the WD. From the location of Quv in a Bq 25 — R 625 
versus Ha — 7^625 color-color diagram and synthetic colors 
of simulated spectra, we estimate that the equivalent width 
of such an emission line is about —70 A. However, the in¬ 
clination inferred from the mass function (Table [1]) and our 
estimate for Mwd (Table [3|) is high (63° to 90°), so eclipses 
are expected if a large cloud of ionized material would be 
present. The absence of radio eclipses in 47 Tuc Q then ar¬ 
gues against the presence of such an extended Ha emission 
region. We note that, given the proximity of the two brighter 
stars, the photometry for Quv should be considered with 
some caution; in particular, the PSFs of relatively bright 
stars in the F625W and Ha images shows Airy disk peaks 
around the central region of the PSF, which coincide with 
the location of Quv and could affect its colors if not properly 
modelled. 

The location of Quv with respect to the He WD cooling 


tracks suggests a mass between ~0.175 and 0.20 Mq. This 
is consistent with the lower limit of the companion mass 
derived from the pulsar mass function (0.18 Mq for an as¬ 
sumed neutron star mass of 1.35 Mq). The derived cooling 
age is very uncertain, as Quv lies in between the tracks for 
the slowly cooling 0.175 Mq models and the rapidly cooling 
0.20 Mq He WDs; these respective masses predict an age 
range of 5.2 — 0.3 Gyr. In addition, the proto-WD phase 
can add up to 1 Gyr to obtain the time that passed since 
Roche lobe detachment. However, a comparison of the WD 
age with the characteristic pulsar spin-down age Tc is diffi¬ 
cult. The observed period derivatives for the 47 Tuc MSPs 
are dominated by the effect of the acceleration of the pul¬ 
sars along the line of sight in the gravitational potential 
of the cluster. As a result, the intrinsic period derivatives, 
and therefore the pulsar characteristic ages, are poorly con¬ 
strained. However, better estimations of the intrinsic period 
derivatives, for stable MSP-WD systems, can be made by 
using the time derivative of the observed orbital periods, 
which essentially measures the acceleration of each binary 
MSP along the line of sight. Details about this technique 
will be given in Freire et al. (in preparation). By using this 
method we obtain Tc > 1.43 Gyr for Quv- 


4.2 47 Tuc S 

Suv, together with Tuv, is the faintest MSP companion 
that we found. This object is clearly visible in our F435W 
images, but in the F625W and Ha images a reliable detec¬ 
tion cannot be made. This is partly the result of a relatively 
bright diffraction spike of a brighter star that passes within 
a few pixels of the center of the PSF of Suv- From com¬ 
parison with the He WD cooling tracks, we derive a mass 
of about 0.2 — 0.25 Mq and an age of about 0.4 — 0.1 Gyr. 
Taking into account its time spent as a ” proto-WD” (up to 
0.4 Gyr), this age range is roughly consistent with the lower 
limit on the pulsar characteristic age. 

47 Tuc S and the isolated pulsar 47 Tuc F are only 0('74 
apart. Both are associated with a single Chandra source, 
W 77, which is likely a blend of (at least) these two sources. 
This explains the relatively large offset between Suv and 
W77 in Tabled 
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4.3 47TucT 

We derive a mass for Tuv of 0.2 — 0.3 M©, in agreement 
with the lower limit of the companion mass (Table [1]). The 
approximate total WD age is about 0.8 —0.1 Gyr. This value 
is roughly consistent with the lower limit on Tc of this MSP. 

Tuv lies in between the CO WD cooling sequence and 
the SMC main sequence in our NUV CMD. The blue color 
of Tuv is apparent from the finding charts in Figure [1] In 
the F390W image Tuv has a close neighbor at a separation 
of only O'.TT. In the F300X image the neighbor, as well as 
stars of similar F390W magnitude, are much fainter. Tuv is 
clearly detected in our F435W images. While also detected 
in F625W and F658N, it is very faint at these redder wave¬ 
lengths; the close neighbor and the contamination of its PSF 
by a weak diffraction spike of another neighboring star, leads 
to relatively large errors in its optical colors (Figure[3]). Nev¬ 
ertheless, the optical counterpart is also clearly blue, and its 
location on the red side of the Ha — i ?625 main sequence 
points at the presence of a strong Ha absorption line. Tuv 
was already identified as the likely counterpart to 47 Tuc T 
by E03a. They detected Tuv in their U images but not in 
V and /, and could therefore only derive an upper limit to 
the color. 


4.4 47 Tuc Y 

Our photometry confirms the blue color of Yuv reported by 
E03a. We estimate that the mass of Yuv is between 0.16- 
0.175 M©, in agreement with the lower limit from the radio 
mass function (Table [T]). The total age estimate, about 2 
Gyr for the cooling time plus 1.3 — 1.9 Gyr for the proto- 
WD phase, makes Yuv, possibly together with Quv, the 
oldest of the five WD companions discussed in this paper. 

Yuv is clearly detected in our E435W image. The PSE 
of the neighboring bright object, which can be seen a little 
offset from the center in the finding charts of Eigure [1] re¬ 
duces the sensitivity for very faint objects. As a result, Yuv 
is not detected in the Re25 and Ha images. 


4.5 47 Tuc U 

EOI found that Uuv is a variable blue object, with a semi¬ 
amplitude variation of 0.004 mag and with a period that is 
consistent with the orbital period derived from radio timing 
(about 0.43 days). Based on a comparison between the U, 
U, an d I photometry for this star and the ISerenelli et al l 
(l200jl models for He WD evolution, its mass and age were 
estimated to be ~0.I7 M© and ~0.6 Gyr, respectively. EOI 
also identified this object as the counterpart of the X-ray 
source WII. EOI indeed found that the He WD companion 
is well inside its Roche lobe, which also indicates there is no 
ongoing accretion that contributes to the X-rays. 

Uuv is the brightest of the five counterparts discussed 
in this work. It is an isolated object, so its UVIS and WEG 
photometry is robust. Uuv is located above and to the left of 
the SMG main sequence in the NUV and optical GMDs, in 
agreement with the results of EOI and E03a. The location 
of Uuv to the red side of the main sequence in the R 625 
versus Ha — Rq 25 GMD points at the presence of a strong 
Ha absorption line, and is indeed in agreement with the He 
WD cooling tracks. 



Figure 4. Finding chart of 47 Tuc I in the NUV. The solid red 
circles are the 3 a match circles for the radio position of the pulsar. 
They are centered on 47 Tuc I. The dashed cyan circles are the 
3 fj match circles of the Chandra source. The dashed white circles 
indicate the blue object that lies at 07l2 from the pulsar. The 
images are l70 x l70 in size. 


The He WD atmosphere models used in this paper are 
computed for a lowe r metallicity and make use of updated 
atmosphere models (|Rohrmann et al.ll20Ilh with respect to 
the ones used in EOI. Erom the updated models and our 
UVIS photometry, we constrain the mass for Uuv to be 
0.16 — 0.175 M©, Teff ~ 11,000 K, and log^ ~ 5.6. This 
is consistent with the results from our WEG photometry 
and with the parameters derived by EOI. Our cooling age 
estimate (0.2 —0.3 Gyr) is somewhat lower than found in the 
latter study (~0.6 Gyr). If we take into account Atproto, the 
total age becomes 2.I-I.6 Gyr; this is still consistent with 
the lower limit to Tc from Table [1] 


4.6 47 Tuc I 

The blue star in the 3 a error circle of 47 Tuc I is very faint. 
The errors in the photometry are considerable, and we plot 
the results from both DAOPHOT and KS2 in Eigure [2] (de¬ 
noted as ol and oIks2, respectively); it is too faint to be 
detected in any of our optical images. This object is too 
blue to be a He WD. While the colors are more consistent 
with a GO WD than a He WD, the former option is not 
likely if the star is truly associated with 47 Tuc I. MSPs 
with GO WD companions typically have longer spin periods 

S ared to systems with He WD companions (|Tauris et al.l 
note the exception of PSR JI6I4-2230). The spin pe¬ 
riod of 47 Tuc I, on the other hand, is very short {Ps = 3.2 
ms; Ereire et al. 2003). We also point out that the mass func¬ 
tion suggests that the companion is a very low-mass object, 
with rric > 0.013 M©. If the blue object is the ~ 0.5 M© GO 
WD companion to 47 Tuc I, this would imply an inclination 
of z ~ 1.8°. None of the aforementioned arguments are con¬ 
clusive grounds for rejecting or confirming the association. 


5 DISCUSSION 

The five binary MSPs in 47 Tuc that we discuss in this work 
are characterized by spin periods below 8 ms, orbital periods 
Pb between 0.43 and 1.2 d, and companions more massive 
than ~0.09 M© (Table [T]). Based on a compilation of avail¬ 
able data for binary pulsars in the Galactic disk, the nature 
of the companion for such binary pulsars is most likely a He 
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WD fsee lTauris et, aLlboid l. We find that our NUV/optical 
photometry for the five MSP counterparts is indeed consis¬ 
tent with the magnitudes and colors predicted by the low - 
metallicity He WD cooling models bv ISerenelli et al] (|2002h . 
In particular, the derived mass constraints lie in the range 
0.16 — 0.30 M©. Our results for the companion of 47Tuc U 
are consistent with those derived by EOl, who compared U, 
y, and / photome try for this star with t he solar-metallicity 
He WD models by ISerenelli et al.l (|200lh . 

The standard formation scenario for an MSP accom¬ 
panied by a He WD involv es a long period of sta ble mass 
transfer. In the models by ISerenelli et al.l (|2002h . the He 
WDs emerge from this phase with a thick H envelope. Cal¬ 
culations for the subsequent evolution of whi te dwarfs pre¬ 
dict a dichotomy in cooling behavioip^ fe.g. [ A lberts et al' 


1996 


2007 


Althaus et al 


Althaus et al 


2001 


; ISerenelli et al.ll20^ : IPanei et al 


20131 ). He WDs with masses ^.0 - 2 M© 


(or equivalently Pb > 1.55 d; see Ivan Kerkwiik et al.|[2005l i 
experience thermonuclear H-shell flashes before they reach 
the cooling track. These flashes leave the WD with a thin H 
envelope, as a result of additional diffusion-induced H-shell 
flashes. After these short-lived episodes, the WD remnant 
cools down on a relatively short timescale. In He WDs less 
massive than ~0.2 M© H-shell flashes do not occur, and 
therefore a thick H envelope is retained. Residual H burning 
can keep up a high temperature over long timescales, result¬ 
ing in much longer cooling times. Figure [2] implies that Suv 
and Tuv cool fast, Uuv and Yuv cool slowly, while for Quv 
the time scale is ambiguous. Unfortunately, the spin-down 
ages of the pulsars in 47Tuc are too uncertain to give an 
independent constraint on the WD age, and hence, mass. 

There are only two binary MSPs in our field whose mass 
functions suggest He WD-like companions but for which we 
have not found a plausible counterpart, viz. 47Tuc E and 
H. Their orbital periods are the longest of the binary MSPs 
in 47Tuc included in our field (~2.3 d). The relation be- 
tween orbital period an d companion mass (see for example 
iTauris Sz Savoniig|l999l i then suggests that the companions 
in these systems are more massive than those in the other 
MSPs. Since a higher mass implies a shorter cooling time, 
the non-detection of the companions to 47 Tuc E and H can 
be understood if they were formed more than ~0.7 Gyr ago 
so that they have dimmed below the detection limit. 

In Figure [5] we show the relation between the mass 
of the He WD and the o rbital period base d on b inary- 
evolution calculations by iTauris Savoni~i3 (|l999l : valid 
for He WDs masses betwe en 0.18 and 0.45 M©) and 
IPe Vito Benvenutol (|2010l i. valid for Pb > 0.25 d. We 
find that the orbital periods and mass estimates for Suv 
and Tuv are mor e consi stent with the relation derived by 
iTauris Savonii j (|l999l i. In contrast, Yuv and Uuv are 
offset from that relation, but their estimated masses lie 
below the mass range to which the relation applies. How¬ 
ever, we observe that Yuv and Uuv are in g ood a greement 
with the relation by IPe Vito Benv enutol (l2010li. Quv is 
consistent with both relations. IPe Vito Benvenutol |2010l i 


The evolutionary models bv llstrate et al.l (l2014l i do not show 
this dichotomy in cooling behavior depending on the occurrence 
of H shell flashes or not, but rather indicate that the cooling 
timescale mainly depends on the mass of the proto-WD. 



He WD Moss (M©) 

Figure 5. Mass-orbital period relations for He WDs. I n gray 
is the relation as calculated by ITauris &: Savoni~i3 (|l999l i. with 
metallicities in the range Z = 0.001 — 0.02 (from right to left) 
and for W D masses in the range 0.18 — 0.45 M©. In orange the 
relation bv IPe Vito fc Benvenutol (|2010li for solar metallicity and 
for orbital periods > 0.25 d. The width of the relation corresponds 
to the uncertainties in its parameters. Overplotted are our results 
for the companions of 47 Tuc Q, S, T, U and Y. The error bars 
denote the possible ma sses for each ob j ect as determined from the 
WD cooling tracks bv ISerenelli et aP (|2002l b We also plot other 
pulsar companions for which Icr uncertainties are less than 10% 
of the WD mass (see ITauris &: van den Henv^ l2014l and refer¬ 
ences therein). The dashed vertical line denote s the lower mass 
limit (0.18 M©) for which the lTanris Sz Savoniid (Il999l) relation is 
valid. The dashed horizont al line shows the minimum or bital pe¬ 
riod (0.25 d) for which the IdI Vito fc Benvennt^ (l2010l) relation 
is valid. 


made different assumptions on the mass transfer process 
and they obtain, for the same or bital periods, lower masse s 
for the WD companions than ITauris Savoniid (| 19991 ). 
They also find that their mass - orbital period relation is 
not sensitive to the initial mass of the accreting neutron 
star. Studies about the assumption on how conservative the 
mass transfer is, have f ound that the dispersion around the 
IPe Vito Benvenutol (|201Q|) relation decreases as the He 
WD mass increases ( De Vito <V Benvenutol [201^ ). 

The five binary MSPs discussed here have X-ray lum i- 
nosities Lx = 2.9—8.9x 10^° erg (|Bogdanov et al.ll2006l ). 
Like for the majority of the 47 Tuc MSPs, it was found that 
their X-ray emission can be described by a thermal (black- 
body or neutron-star hydrogen atmosphere) spectrum as ex¬ 
pected for emission from the heated magnetic polar caps of 
the neutron stars. Using the L, Te// and mass estimates 
from Table [3l and the orbital periods from Table [U we have 
calculated the WD radii (Ru^d) and (assuming a neutron- 
star mass of 1.35 M©) the corresp onding radii of t he Roche 
Lobe {Rl,wd] using the formula of |Paczvhskilll97ll ) . We in¬ 
deed find that all five WD companions fit well within their 
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Roche lobe {Rl,wd/Rwd ^ 5) as already found for 47Tuc 
U by EOl, so we indeed expect that accretion does not con¬ 
tribute to the X-ray luminosity. 

For the five systems we studied, we calculate the contri¬ 
bution of the MSP radiation to the observed NUV luminos¬ 
ity of the WD. We used the value s of the spin-down lum i- 
nosity (E) for each MSP given bv iBogdanov et aD (|2006h . 
Assuming isotropic emission from the neutron star we cal¬ 
culate the energy flux intercepted by the WD (Ei) using our 
estimates of Rwd and the separation between the stars. For 
Quv, Suv and Tuv we obtain log Ei (Lq) < —4 which is well 
below their NUV luminosities. This result indicates that ir¬ 
radiation is not driving the luminosities of these systems. On 
the other hand, for Yuv and Uuv we obtain —4.2 <\og Ei 
{Lq)< —2.3 which is a closer value to their respective L, 
implying that heating might be present, though not domi¬ 
nating the NUV luminosity. This result can be understood if 
we take into account that these two companions have larger 
Rwd and are closer to their MSPs than the other three 
WDs. 

We note that all the companions have cooling ages < 6 
Gyr (including age uncertainties), which are considerably 
less t h an the age o f the c luster (9.9 ±0.7 Gvr. lHansen et aD 
I 2 OI 3 I I . iRasio et al.l (|2000l l proposed that the formation of re¬ 
cycled MSP systems happened mostly during the early life 
of the cluster. However, our results may suggest that dy¬ 
namical interactions during the later evolution of the clus¬ 
ter have a more important role in the formation of these 
systems (for an evoluti onary scenario in this context see for 
example lFregeaull2008h . 

The orbital periods of the five binary MSPs studied in 
this paper are all longer than the time span of the NUV ob¬ 
servations (~9 h); light curves constructed from these data 
alone give incomplete phase coverage. We defer a variability 
study, based on these and other HST data sets, to a future 
paper. 
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